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"  The  FACT  (Fast  Asymptotic  Coherent  Transmission)  Model  is 
the  new  Navy  Interim  Standard  Transmission  Model  for  Ocean 
environments  which  may  be  treated  within  the  context  of  a  single 
sound-speed  profile  and  a  flat  bottom.  It  is  a  ray-acovntica 
model  designed  for  the  computation  of  transmission  loss  as  a 
function  of  range  and  frequency  at  fixed  depth.  The  classical 
ray  treatment  has  been  augmented  with  higher  order  asymptotic  ^ 
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4.  PROGRAM  STRUCTURE  AND  FLOW 


This  section  is  divided  into  four  subsections: 

4.1  Description  of  inputs 

4.2  Sequence  of  calculations 

4.J  Description  o i  outputs 

4.4  Detailed  flow  charts  and  descriptions  of  each 
of  the  basic  programs  and  subroutines  in  the 
FACT  package. 

These  sections  were  taken  largely  from  the  Ocean  Data 
Systems,  Inc.  final  report  on  FACT  model  development  for 
ai:sd.  • 

The  basic  flow  of  the  FACT  model,  exclusive  of  input 
and  output  is  summarized  in  Figure  4-1,  and  is  considerably 
expanded  upon  in  sections  4.2  and  4.3.  The  complete  pro¬ 
gram  listing  and  punched-card  deck  are  available  from  AESD 
upon  request. 

The  objective  of  FACT  is  to  estimate,  by  using  ray¬ 
tracing  techniques,  the  acoustic  transmission  loss  in  a 
single-profile,  flat-bottom  ocean  environment,  as  a  function 
of  range  and  frequency.  Additionally,  if  requested,  FACT 
will  produce  the  arrival  angles  (at  the  receiver)  of 
individual  ray  paths,  again  as  a  function  of  range.  Trans¬ 
mission  loss  (dB  re  1  yard)  is  tabulated  in  a  single  array 

•c.  L.  Baker,  FACT  Transmission  Loss  Model  Development  for 
Acoustic  Environmental  Support  betachment  -  ^lnal  Report, 
Ocean  Data  Systems,  Inc.,  June  30,  1973  on  Contract 
N00014-73-C-0131. 
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of  dimension  250x6  at  up  to  250  equally  spaced  range  points 
for  each  of  one  to  six  frequencies.  Arrival  information  is 
written  to  an  auxiliary  (tape  or  disk)  file  as  individual 
records  containing  fields  for  range,  angle,  and  inten¬ 
sities  at  up  to  six  frequencies. 

As  indicated  in  the  documentation  inc’uded  as  part  of 
the  FACT  Handout,  the  primary  component  of  the  FACT  Package 
is  a  single  subroutine  FACTTL,  which  may  be  incorporated 
into  any  of  a  number  of  complete  programs  requiring  an 
estimation  of  transmission  loss  versus  range  and  frequency. 
One  example  of  a  stand-alone  program  is  included:  TLOSS , 
a  program  which  reads  input  parameters  from  cards,  calls 
on  FACTTL  for  losses,  and  prints  and/or  plots  the  results. 
This  program  is  primarily  useful  to  analysts  requiring  a 
small  number  of  runs  as  part  of  a  design  program  on  a 
demand  basis. 

Two  additional  transmission  loss  models  may  be  used  to 
supplement  FACT  in  those  cases  where  a  full  FACT  solution 
is  liable  to  result  in  excessive  running  times.  These 
models,  SHALTL  and  HFCHTL,  are  designed  specifically  to 
approximate  the  results  of  a  complete  FACT  solution  in 
shallow-water  transmission  and  half-channel  transmission 
respectively.  Subroutine  HFCHTL  is  an  integral  component 
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of  FACTTL  in  that  the  output  of  HFCHTL  is  supplemented  by 
the  output  of  FACTTL  for  the  direct  and  bottom-reflected 
paths.  On  the  other  hand,  in  order  to  employ  subroutine 
SHALTL  a  modification  to  TLOSS  is  required  (e.g. ,  replacing 
the  call  to  FACTTL  by  a  call  to  SHALTL).  Care  should  be 
exercised  in  using  both  of  these  models  as  both  serve 
only  as  quick-running  alternatives  to  the  normal  FACT  pro¬ 
cessing.  Additionally,  the  HFCHTL  model  requires  further 
care  m  use  m  that  it  is  valid  only  for  the  specific 
frequencies  and  source/receiver  combinations  contained 
within  the  listing. 

The  documentation  which  follows  is  designed  to  supple¬ 
ment  the  user-oriented  material  presented  in  the  FACT 
Handout.  It  is  intended  for  computer  programmers,  and  is 
primarily  concerned  with  the  overall  structure  and  inter¬ 
relation  of  the  various  components  of  the  FACT  Package. 

The  organization  of  this  material  follows  that  of  the  FACT 
programs  components  is  presented  in  Section  4.4. 

The  FACT  model  itself,  written  entirely  in  FORTRAN  IV, 
is  invoked  by  a  single  call  to  subroutine  FACTTL.  Core 
requirement,  excluding  input  and  output,  but  including 
all  other  FACT  and  system  computational  routines,  is 
approximately  8,400  decimal  ^(20,300  octal)  cells  on  the 
CDC  6400/6600. 


4-3 


f 


4. 1  Inputs 

The  inputs  to  FACT  are  primarily  geometrical  and 
environmental  in  nature.  They  are: 

.  A  sound  velocity  profile:  the  speed  of  sound  is 
specified  at  each  of  up  to  50  points  from  the 
surface  to  the  bottom.  Depth/velocity  pairs 
are  in  feet  and  feet  per  second  or  meters 
and  meters  per  second. 

.  Surface  and  bottom  conditions:  wave  height  in 
feet;  FNWC  bottom  class. 

.  Source  and  receiver  positions:  depths  in  feet. 

.  Frequency  information:  frequencies  in  Hertz; 
coherency  flags. 

Range  information:  number  and  spacing  of  range 
points  in  nautical  miles. 

4 . 2  Sequence  of  Calculation 

In  the  following  presentation,  only  the  most  signifi¬ 
cant  steps  in  determining  transmission  loss  are  outlined; 
many  computational  steps,  such  as  the  calculation  of  con¬ 
stants  and  other  factors  essential  to  the  calculation  are 
covered  in  detail  in  the  sections  dealing  with  individual 
subroutines.  Some  liberties  have  been  taken  in  describing 
the  sequences  of  calculations,  but  it  is  essentially  as 
follows : 
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Profile  correction:  The  profile  points  are  corrected 
to  take  account  of  spherical  earth  geometry.  Note 
that  depths  and  depth  indices  are  often  used  inter¬ 
changeably  (as  appropriate)  in  discussing  FACT. 

Axis  location:  The  deep  sound  channel  axis,  if  any, 
of  the  profile  is  located,  and,  under  certain 
conditions  the  source  and  receiver  depths  are 
altered  to  allow  simulation  of  axis-to-cxis  trans¬ 
mission. 

Profile  augmentation:  The  source  and  receiver  depths 
are  inserted  in  the  profile  as  explicit  points, 
altered  slightly,  if  necessary,  to  avoid  equal 
velocities  at  the  two  depths. 

Source  selection:  For  the  remainder  of  the  FACT 

calculat ions ,  the  source  and  receiver  depths  are 
selected  from  the  two  deptM  of  interest  by 
assigning  the  depth  with  the  leaser  sound  velocity 
as  the  source  depth  --  the  starting  point  for  the 
ray  tracing  process.  The  index  of  this  point  in 
the  profile,  Ki,  is  used  heroarter  to  designate 
tne  source  (or  source  depth):  the  other  ii»dex,  K2 , 
is  used  to  designate  the  receiver  or  receiver  depth. 
(Cars  Is  taken  to  distinguish  the  "true"  source 
from  the  ray-tracing  source  where  the  distinction 
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ia  important.) 


.  Geometry  factors:  A  number  of  flags  are  set  (at 
various  points  throughout  the  program)  to  in¬ 
dicate  various  geometrical  relationships  be¬ 
tween  source  and  receiver. 

.  Low-f requency  effects:  The  WK^  phase  factors  for 
low  frequency  cutoff  are  calculated. 

.  Ray  selection:  The  angles  of  the  rays  to  be  traced 
from  K 1  to  K2  are  selected,  grouped  into  one  or 
more  families.  The  selection  is  based  on  the 
velocity  profile,  source,  and  receiver  depths. 

Rays  are  chosen  so  that  within  each  family,  an 
analytical  fit  of  Range  versus  Angle  can  be  made, 
thus  smoothing  and  retaining  legitimate  caustics 
while  removing  false  caustics;  the  functional 
form  of  the  fit  will  vary  with  family  type.  If 
the  profile  and  associated  source  and  receiver 
depths  lead  either  to  more  than  20  families  or  100 
rays,  processing  is  terminated  and  the  transmission 
loss  array  is  returned  with  zero  values  for  all 
entries . 

.  Ray  tracing  notes:  Because  the  environment  is  single 
profile,  flat  bottom,  any  ray  which  is  traced 
exhibits  a  periodicity  over  the  range  of  interest 
and  is  actually  traced  for  only  a  single  such  cycle. 

#  Wentzei-Kramers-Brillouin  approximation  -  see  L.  M. 

Brekhovskikh,  Waves  in  Layered  Media,  Academic  Press,  New 

York,  1960b,  pp  452-3. 
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Within  a  single  period,  each  ray  angle  selected 
may  represent  one,  two,  or  four  paths  (arrivals) 
from  source  to  receiver,  per  period,  depending 
upon  the  geometry  involved:  the  source  angle  may 
be  positive,  negative,  c  r  L  -th,  and  ray  either  1» 
reflected  cr  refracted  at  :a,<  j\*«  r,  ri 

cross  and  re-cross  the  receiver  depths.  Cac). 
arrival  is  assigned  a  number,  termed  the  arrival 
order,  indicating  the  number  of  deep  cycles  which 
the  ray  has  undergone;  arrival  order  corresponds 
to  the  direct  path.  The  ordering  of  arrival  ranges 
within  a  single  period  or  cycle  of  a  ray,  and  the 
assignment  of  these  ranges  (plus  multiples  of  the 
period)  to  individual  arrival  orders  is  determined 
and  controlled  by  a  number  of  flags  and  indices 
which  are  detailed  elsewhere. 

Path  combinations!  Depending  upon  the  geometries  in¬ 
volved,  two  or  four  paths  from  the  source  to  the 
receiver  may  be  combined  into  a  single  path  of 
doubled  or  quadrupled  Intensity.  This  may  happen 
for  instance,  if  the  source  is  so  close  to  the  sur¬ 
face  that  the  downgoing  ray  and  the  surface-reflected 
upgoing  ray  are  essentially  parallel  and  arrive  at 
the  receiver  at  essentially  identical  angles.  A 
similar  geometry  may  apply  at  the  receiver,  or  at 
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both  the  source  and  receiver.  The  processing 
of  such  combined  rays  is  controlled  by  several 
flags  indicating  the  simplified  situation. 

Family  processing:  Each  family,  in  turn,  is  processed 
to  determine  its  contributions  to  the  total  in¬ 
tensity  arising  from  the  source,  and  equivalently, 
overall  transmission  loss.  At  the  same  time, 
individual  arrival  angles  and  intensities  are 
written  to  a  separate  file  for  later  processing. 
This  loop  on  families  is  the  major  processing 
loop  in  FACT-  Within  each  family,  the  atero- 
order  or  direct  path  is  processed  first.  This 
is  followed  by  a  loop  in  which  subsequent  orders 
1,  2,  3,...  are  processed  in  turn.  There  is  no 
fixed  limit  on  the  number  of  orders  which  FACT  may 
be  required  to  process;  this  loop  is  terminated 
only  when  significant  intensity  no  longer  is  being 
contributed  by  any  path  within  the  family  at  the 
ranges  of  interest. 

Half-channel  note:  When  a  half-channel  case  has 

been  flagged  on  input,  only  the  direct  and  bottom 
and  surface-reflected  arrivals  are  processed  as 
outlined  above.  In  these  cases,  the  non-direct 
path,  non-bottom  and  surface  reflected  contribution 
to  intensities  are  approximated  and  added  by  a 
separate  half-channel  model. 
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.  Final  procotingt  When  ell  families  have  been  pro¬ 
cessed,  surface-duct  contributions  if  present, 
are  added  to  those  intensities  already  calculated, 
and  then  converted  to  transmission  loasea  (re  one 
yard) . 

Processing  of  an  arrival  order  of  a  family  of  raye  consists 
of  the  following  steps: 

The  arrival  ranges  for  each  of  the  (one  to  four)  paths 
with  this  order  are  calculated. 

For  each  "ath.  the  coefficients  and  parameters  re¬ 
quired  to  express  range  as  a  function  of  ray  angle 
are  calculated.  Any  one  of  four  possible  functional 
forms  is  used  according  to  family  characteristics. 

It  the  range  Intervals  for  all  four  path^  exceed  the 
maximum  range  of  interest,  processing  of  arrival 
orders  for  the  family  is  terminated. 

.  Subroutine  XNSTOK  (or  CUSP  If  applicable)  is  called 
to  calculate  and  add  the  intensity  arising  from 
each  (smoothed  and  fitted)  path  to  the  transmission 
loss  array  at  each  rang#  point  for  such  frequency. 

.  Zf  the  intensities  from  all  four  patha  drop  balow  a 
specified  minimum  value,  procsssing  of  arrival 
orders  for  the  family  la  terminated. 


Processing  of  one  path  of  en  errivel  order  by  ZtiSTOW  or  CUSP 
consists  of  the  following  steps » 


The  type  of  fit  of  range  versus  ray  angle  is  examined 
to  determine  whether  or  not  a  caustic  exists,  and 
to  find  the  minimum  and  maximum  ranges  at  which 
contributions  to  total  intensity  are  made. 

If  this  range  interval  is  bevond  the  range  of  interest, 
processing  of  the  path  is  terrir.ated . 

At  each  applicable  range  point,  the  number  of  arrivals 
(rays)  is  calculated:  zero  indicates  the  shadow 
zone  of  a  caustic,  one  or  two  indicates  an  illumin¬ 
ated  region. 

The  intensity  contribution  from  each  ray  is  added  to 
the  transmission  loss  array  for  each  frequency  at 
the  range  being  processed.  The  intensity  is  com¬ 
puted  as  an  analytic  function  of  range,  frequency, 
and  the  values  of  ray  anqlc  and  the  derivations  of 
range  with  ray  angle  at  this  range;  the  latter  are 
obtained  by  examination  of  the  range  versus  ray 
anq 1 e  fit. 

The  calculated  intensities  are  modified,  if  required  by 
factors  reflecting  coherent,  semi -coherent ,  or  inco¬ 
herent  path  addition,  shadow-zone  fall-off,  low- 
frequency  cutoff  effects,  and  bottom-bounce  losses 
as  applicable. 

If  flagged,  range,  arrival  angle,  and  intensity  infor¬ 
mation  is  written  to  an  external  file. 


Whan  all  range  point*  have  been  processed*  a  flag  Is 
set  to  indicate  if  the  minimum  range  of  the  path 
has  exceeded  the  range  of  Interest,  or  If  the 
contribution  to  Intensity  has  dropped  below  a 
specified  minimum  value. 


4 . 3  Outputs 

The  primary  output  from  subroutine  FACTTL  is  an  array*  TL, 
of  dimension  250x6,  giving  transmission  loss  (in  dB  re  1  yard) 
at  each  of  the  range  points  and  frequencies  specified  as  In¬ 
put  parameters.  If  the  ray  selection  process  results  either 
in  too  many  rays  (>  100)  or  too  many  families  (>  20) ,  a  two- 
line  message  will  be  printed  Indicating  this  condition,  and 
the  TL  array  will  contain  vc.oa  .>•  the  specified  ranges  and 
frequencies . 

If,  in  addition,  arrival  information  is  desired,  the  < - 

parameter  IAR  should  be  set  indicating  the  FORTRAN  unit  to  wr.i  h 
arrival  records  should  be  written.  The  forsuit  of  these  re¬ 
cords  is  detailed  in  the  description  of  subroutine  INSTOR. 

The  output  file  is  never  positioned  nor  are  any  file  marxs 
written;  these  functions  are  delegated  to  the  calling  program. 

Optional  debugging  output  is  also  available  via  an  input 
parameter;  care  should  be  taken  in  setting  this  flag,  as  the 
potential  exists  for  many  hundreds  of  pages  of  output  from  a 
single  run. 
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4 . 4  Components  of  the  FACT  Package 

This  section  consists  of  a  description  of  each  of  the 
components  of  the  FACT  Package  --  the  driver  programs,  the 
main  computational  routine,  and  auxiliary  subroutines  and 
functions.  For  each  of  these,  the  following  material  is 
included : 

.  A  brief  description  of  the  function  of  the  component 
in  the  Model. 

.  Equations  used  by  the  component  when  these  are  not 
immediately  evident  from  the  function  of  the 
component  and/or  the  program  listing. 

.  Parametric  and  common  input  and  output  variables. 

.  Flow  charts,  expressed  in  physical  terms  to  the 

greatest  extent  possible,  for  the  major  programs 
and  routines  of  the  Model. 

.  Additional  material,  as  applicable,  to  present  the 

details  of  the  program  logic  not  included  m  the 
flow  charts  for  the  routine. 


INSERT  SOURCE  AND  RECEIVER  IN  PROFILE 


DETERMINE  RAYS  TO  BE  TRACED 


LOOP  ON  RAY  FAMILY 


TRACE  RAYS  FOR  FIRST  CYCLE 


PROCESS  DIRECT  PATHS  -  COMPUTE  NPATHS 


LOOP  ON  NPATHS 


FIT  R {0 ) 


PROCESS  TL  (R)  IN  INSTOR 


END  LOOP  ON  NPATHS  fc  PROCESSING  OF  DIRECT  PATHS 


PROCESS  ALL  SUBSEQUENT  ARRIVAL  ORDERS  -  RECOMPUTE  NPATHS 


INCREMENT  ALL  ARRIVALS  BY  ONE  PERIOD 


LOOP  ON  ..PATHS 


DOES  FAMILY  CONTAIN  A  CUSP ED  CAUSTIC 


ES 


- . -,->--w.^y« 


^.ir^iJi*KSt26tsm£iai31Vti}r.!r<f*  ...... 


PROGRAM  TLOSS 


TLOSS  is  the  driver  which  reads  card  inputs,  invokes 
FACTL  to  compute  transmission  loss  (and  arrival  structure)  versus 
range  and  frequency,  and  prints  and/or  plots  the  results. 

CARD  INPUTS 

The  card  input  formats  and  sequence  are  detailed  in  the 
FACT  Handout,  pages  5-9  throuqh  12.  See  also  subroutine  RDPROF, 
page  4-17. 

PRINTED  OUTPUTS 

The  formats  for  tabulated  and  plotted  printer  outputs  are 
detailed  in  the  FACT  Handout,  pages  5-14  through  33.  See  also 
subroutine  TLHEAD,  and  PLOTTL,  pages  4-19  and  20.  TLOSS  also 
prints  the  data  on  the  input  cards  as  they  are  read. 

COMMON  USAGE 
None 


NOTE 

A  complete  FORTRAN  listing  of  TLOSS  is  included  in  the 
portion  of  the  FACT  Handout  included  in  this  report,  pages  5-43 
through  47. 
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Continue  TLOSS 


7S 

Call  FACTTL  to  calcu¬ 
late  transmission  loss, 
and  (if  requested) 
arrival  structure 

r 

1  Yes 

Was  case  aborted  by 

FACTL? 

Print  losses  vs.  range 
if  tabulated  listing 
has  been  requested 

f 


Call  TLHEAD  and  PLOTTL 
if  line-printer  plots 
of  losses  vs.  frequency! 
have  been  requested 

f— 


No 

Was  arrival  structure 

calculated? 

l!1 

Road  arrival  structure 
file,  store  in  trans¬ 
mission  loss  array  as 
integer  angle  less  than] 
25*  vs.  range 


res 


Call  TLHEAD  and  PLOTTL 
to  plot  arrival  angles 
vs.  range  on  line 
printer _ 
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TLOSS  2 


SUBROUTINE  RDPROF 

RDPROF  is  called  from  TLOSS  to  reed  profile  data  ."rom 
carda,  and  convert,  if  necessary,  to  FACT  units.  Card  formats 
are  detailed  in  the  FACT  Handout.  Input  data  is  printed  as 
it  is  read;  see  pages  5-14,  5-23and  24,  and  5-33. 

PARAMETER  INPUT 

N  t  No.  of  Profile  Points 


CARD  INPUTS 

A)  N  Positive: 

Z  Profile  Depths  (4/'Card)  ,  meters  or  feet 
C  Profile  Velocities  (4/Card),  mtra/sec  or  ft/sec 

B)  N  Negative: 

D  Profile  Depths,  r?t  »rs 

T  Profile  Temf  natures ,  degrees  Centigrade 

S  Profile  Salinities,  parts/thousand 


PRINTED  OUTPUTS 

The  data  read  from  cards  is  printed  along  with  the  results 
of  any  conversions  performed.  One  of  three  applicable  formats 
is  used: 


1)  Depths  in  feet  and  velocities  in  feet  per  second. 

2)  Depths  in  meters  and  velocities  in  meters  per  second, 
plus  depths  in  feet  and  velocities  in  feet  per  second. 

3)  Depths  in  meters,  temperatures  in  degrees  Centigrade, 


& 
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RDPRQP  (Cont|d) 


salinities  in  parts  per  thousand,  plus  depths  and 
and  velocities  as  in  2). 

PARAMETER  OUTPUTS 

Z  Array  {50)  of  Profile  Depths,  feet 

C  Array  (50)  of  Profile  Velocities,  ft/sec 

COMMON  STORAGE 
None 
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SUBROUTINE  TLHEAD 

TLHEAD  it  called  from  TLOSS  for  the  purpose  of  printing 
any  of  three  different  heading  formats;  these  are  detailed  in 
the  FACT  Handout. 

PARAMETER  INPUTS 
TITLE 
S 
R 
IF 
JC 
IX 
NF 

ITYP 

PRINTED  OUTPUTS 

A)  ITYP  -  0: 

Page  Heading  for  Tabular  Transmission  Loss  vs.  Range 

B)  ITYP  -  1: 

Page  Heading  for  Plot  of  Transmission  Less  vs.  Range 

C)  ITYP  -  2: 

Page  Heading  for  Plot  of  Arrival  Angle  vs.  Range 

CQMH6R  storage 
Hone 


80  Character  Title  (8  words,  10  thar/vord) 

Source  Depth,  feet 

Receiver  Depth,  feet 

Array  (6)  of  Frequencies,  Hz 

Array  <61  of  Coherency  Flags,  (0,1  or  2) 

Array  (6)  of  Characters  to  be  Used  for  Plottinq 
No.  of  Frequencies 
Heading  Type,  0,1,2 


♦-It 


SUBROUTINE  PLOTTL 


PLOTTL  is  called  by  FACTTL  for  the  purpose  of  plotting, 
using  the  line  printer,  either  transmission  loss  versus  range 
or  arrival  angle  versus  range.  The  plot  format,  examples  of 
which  appear  in  the  FACT  Handout,  consists  of  an  array  of 
printer  positions  121  columns  wide  and  51  columns  high, 
augmented  by  horizontal  and  vertical  scale  information. 

Input  is  taken  from  an  array,  TL,  of  dimension  250  (range 
points)  by  6  (frequencies  or  angles).  Up  to  250  range  points 
and  up  to  6  frequencies  or  angles  may  be  plotted.  Either  a 
separate  (input)  or  a  single  (default)  plotting  character  may 
be  used  for  each  frequency  or  angle;  this  is  controlled  by  flag 
IX.  Alternate  range  points  only  are  plotted  when  the  number 
of  range  points  exceeds  120.  MNDB  is  the  minimum  value  of  the 
plot:  If  greater  than  zco,  the  abscissa  values  increase  down¬ 
ward  from  MNDB  to  MNDB  ♦  50;  if  less  than  zero,  the  abscissa 
values  decrease  downward  from  ABS(NNDB)  to  ABS(MNDS)  -  50.  The 
units  of  the  incremental  range.  DR,  are  arbitrary,  but  such  that 
the  maximum  range  is  9999.9  inn)  or  loss. 

PARAMETER  INPUTS 

Nft  Number  of  range  points,  £  250 
DR  Incremental  range,  nautical  miles 

NF  Number  of  frequencies,  <6 

IX  Array  (6)  of  plotting  characters;  or  zero 
MNDB  Minimum  value  of  output  plot,  dB  or  degrees 

TL  Array  (250,6)  of  value*  to  be  plotted,  dB  or  degrees 
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PLOTTL  (Cont'd) 

PRINTED  OUTPUT 

Lin*  printer  plot  of  the  value*  in  TL  a*  a  function  of 
range,  with  up  to  aix  plotted  values  at  each  of  up  to  120 
range  points. 


COMMON  STORAGE 
None 


SUBROUTINE  FACTTL 


FACTTL  is  the  subroutine  which  is  to  be  invoked  in 
order  to  calculate  transmission  loss  (and  arrival  angle)  vs. 
range.  All  input  and  output  is  directed  by  the  values  of  the 
parameters  of  the  call  to  FACTTL;  the  common  storage  of  FACT 
and  its  subroutines  must  be  kept  separate  from  that  of  any 
calling  routine.  A  condensed  description  of  the  major  pro¬ 
cessing  done  by  FACTTL  is  found  on  pages  5-5  through  7. 

PARAMETER  INPUTS 

YS  Source  depth,  feet 

YR  Receiver  depth,  feet 

FREQB  Array  (6)  of  frequencies,  Hr 

IC  Array  (6)  of  coherency  flags;  0,1  or  2 

WHf  Wave  height,  feet 

IB  FNWC  bottom  type,  1-9 

NPTSPP  No.  of  points  in  sound  velocity  profile,  <  50 

IL  Index  of  surface  layer  in  profile,  £  50 

YPP  Array  (50)  of  profile  depths,  feet 

CPP  Array  (50)  of  profile  velocities,  feet  per  second 

NR  No.  of  range  points,  4  250 

DR  Incremental  range,  feet 

IARVTP  File  unit  for  arrival  output  if  f  0 

IP  Debugging  print  flag,  0  or  1 
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FACTTL  (Cont'd) 


PARAMETER  OUTPUTS 

TL  Array  (250,6)  of  transmission  loaf  va.  range 

and  frequency,  dB 

FILE  OUTPUT 

Unit  IARVTP  One  record  for  each  arrival  angle  at  each 

range  point,  including  path  loaa  va.  frequency. 
(See  INSTGR  for  format,  page  4-79) 


COMMON  VALUES  CALCULATED 
/RANGEL/  N RANGE 

NFREQ 
IFQNIN 
FREQ 
FREQK 
RANGE 

/FLAGS/  IGTYP 

THMIN 
THMAX 
CONST 
C1C2 

CBC2 

ICOH 

IRSR 

NBOT 

|- 


No.  of  range  pointa  (£250) 

No.  of  frequencies  {£  6) 

Index  of  lowest  frequency  6 
Array  (6)  of  (radian  frequencies) **-1/3 
Array  (6)  of  frequencies  in  KHz 
Arra,.  2  "  of  range  values,  feet 
Tvpe  of  family  being  processed  1-4 
Angle  giving  minimum  range  in  fit.  Rad 
Angle  giving  maximum  range  in  fit,  Rad 
Const* (Receiver  Velocity) **1/3 
Ratio  of  source  depth  velocity/ 
receiver  depth  velocity 
Ratio  of  bottom  velocity/receiver 
depth  velocity 

Flag  to  Indicate  combination  of  arrivals 
Flag  to  indicate  surface  reflection 
No.  of  bottom  reflections 
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FACTTL  (Cont  *d) 


IBTYP 

FNWC  bottom  type 

IFLAC 

Array  (6)  of  coherency  flags 

DK 

Array  (2,6)  of  semi-coherent  phase  factors 

/RAYZER/ 

IREFRZ 

Flag  indicating  grazing  family 

THETMB 

Critical  angle  for  bottom  type,  radians 

THBINC 

Incremental  angle  for  bottom  rays,  radians 

/INIT/ 

AK 

Vertex  velocity  of  ray  to  be  traced 

SINTHO 

Sine  of  initial  angle  of  ray  to  be  traced 

I  ML 

Index  of  mixing  layer  in  input  profile 

/INPUTS/ 

I  RAY 

Number  cf  rays  (in  TH)  in  family  being 

processed  ( < 100) 

KP 

Index  (in  R)  of  arrival  being  processed 

(1-4  ) 

NP 

Index  of  path  being  processed  (1-45 

NORDEH 

Arrival  order  being  processed 

NG 

Index  of  group  (family)  being  processed 

'<20) 

IPAR 

Flag  for  type  of  family  being  processed 

(1  or  2} 

R 

Array  (100,4)  of  arrival  ranges  vs.  angle 

and  arrival,  ft 

TH 

Array  (100)  of  ray  angles  in  family 

being  processed,  rad 

/FITS/ 

THF 

Array  (3)  of  angle*  of  fit  of  R  vs.  • 

A 

Array  (S)  of  coefficient*  of  fit 

4-24 


/CUSPCM/ 

/AUTCOH/ 


XN1N  Ordinate  value  for  min  value  of  fit 

XMAX  Ordinate  value  for  max  value  of  fit 

RMXN  Minimum  value  of  range  in  fit 

RMAX  Maximum  value  of  range  in  fit 

RANKIN  Minimum  range  value  to  be  fit 

CCUSP  Velocity  at  cusp  depth,  ft/sec 

FNMIN  Min  no.  of  range  pti/surf  image  cycle  ■  8/3 

FNMAX  Max  no.  of  range  pts/surf  image  cycle  ■  6 

FNXI  Reciprocal  of  FNMAX-FNMIN 

FNCYC  Array  (2,65  of  cycles  of  phase  difference 


NOTE 

Following  the  flow  charts  for  FACTTL,  additional  pages  are 
presented  which  elaborate  on  the  geometry  of  the  ray  paths 
which  are  used  in  the  FACT  trans  *  i  >n  .  ss  calculation,  and 
the  values  of  flags  and  indK.  .>  required  to  sequence  the 
corresponding  path  ranges  as  a  function  of  family,  arrival 
order  and  path  index. 
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Continue  FACTTL 


Exit  FACTTL 
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FACTTL  2 


FACTTL  -  Loop  on  Families 


A-  •  _  t  C  —  —  X _ J  %  ■ 

AitAc^axA-xc;  a  am  a  Ay  paiajn'” 

eters 


Loop  over  all  rays  in 
family  IR  ■  NR1 ,  NR2 


Call  RANGER  to  calcu 
late  ray  period  and 
source-receiver  range** 
of  first  two  arrivals 
of  ray _ 


_L 

Calculate  source- 
receiver  ranges  of 
third  and  fourth 
j  arrivals  of  ray 


Calculate  minimum  and  maxi¬ 
mum  ray  angles  for  family; 
calculate  surface  reflec- 
ticn  parameters 


For  surface-reflected 
family,  compute  phase 
relationships  of  up-  and 
down-going  rays 

Determine  geometry  for 
pairing  of  arrivals  and 
set  ICOH  flat 
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FACTTL  3 


Continue  FACTTL  -  Loop  on  Families 


FACTTL  -  Loop  on  paths  of  higher  arrival  orders 
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FACTTL  7 
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FACTTL  -  Coherence  Gaometry  for  Combining  Ray  Path* 
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SUBROUTINE  SHALTL 


SHALTL  ia  a  self-contained  model  for  estimating  shallow- 
water  (less  than  1000  feet)  transmission  loss  for  bottom- 
class/f requency  combinations  permitting  perfect  reflection 
of  rays  at  grazing  angles  on  the  bottom  less  than  some 
critical  angles.  A  trivial  modification  to  TLOSS  will  cause 
SHALTL  to  be  called  in  place  of  FACTTL  for  these  combinations 
which  may  bo  followed  for  a  call  to  FACTTL  for  the  remaining 
cases.  Note  that  a  single  frequency  is  processed  by  each  call 
to  SHALTL. 

SHALTL  assumes  a  homogeneous  (uniform  sound  velocity) 
medium  and  includes  an  average  surface  image  interference 
effect.  A  simplified  bottom-loss  approximation  assumes  a 
bottom  which  is  perfectly  reflecting  to  a  critical  angle 
THCDG.  Above  this  angle,  the  bottom  has  a  constant  loss  FL90 
for  the  first  order  path  and  is  perfectly  absorbing  for  higher 
order  paths.  The  constants  THCDG  and  FL90  are  chosen  from 
3x3  arrays  as  a  function  of  bottom  type  (1-3)  and  frequency 


band  (0- 

150  Hz,  151-699  Hz,  and 

700- 

1000  Hz). 

PARAMETER  INPUTS 

YS 

Source  depth,  feet 

D 

Bottom  depth,  feet 

YR 

Receiver  depth,  feet 

NR 

No.  of  range  points,  <250 

F 

Frequency,  Hz 

DR 

Incremental  range,  feet 

ZB 

Bottom  type,  1-3 
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SHALTL  (Cont'd) 


PARAMETER 

TL 


OUTPUT 


Array  (25;)  of  transmission  loss  versus  range, 
dB  re  3  yard 
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SUBROUTINE  INSERT 


INSERT  is  called  by  PACTTL  at  the  beginning  of  each  case 
to  procoss  the  input  sound  velocity  profile.  INSERT  corrects 
all  depths  and  velocities  to  account  for  spherical-earth 
effects,  and  ensures  that  explicit  points  for  source  and 
receiver  depths  are  inserted  in  the  profile  and  that  they  are 
at  points  of  unequal  sound  velocity.  Prior  to  inserting 
these  sourco  and  receiver  depths  in  the  profile,  INSERT  calls 
AXIS  to  account  for  axis-tc-axis  transmission;  thia  call  may 
result  in  the  source  and  receiver  depths  being  changed.  For 
ray-tracing  purposes,  INSERT  chooses,  from  these  two  depths, 
the  point  with  the  lesser  sound  speed  as  the  source  for  ray¬ 
tracing  purposes;  throughout  the  remainder  of  the  FACT  program, 
the  term  source  refers  to  this  point.  Three  indices  are  set 
by  INSERT  as  the  result  of  this  processing:  Kl,the  index  of 
the  ray-tracing  source;  K2,  the  index  of  the  ray-tracing 
receiver;  and  KRC  («1  or  2)  indicating  which  of  these  is  the 
depth  at  which  arrival  angles  are  to  be  calculated. 

INSERT  constructs  the  final  sound  velocity  profile  in 
common  area  /VEbOC/  and  the  gradients  of  this  profile  in  coounon 
area  /GRADS/. 

PARAMETER  INPUTS 

Y1  Source  depth,  feet 

¥2  Receiver  depth,  feet 

NPTS  No.  of  points  in  sound  velocity  profile  50 
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INSERT  (Cont’d) 


YX  Array  (50)  of  profile  depth*,  feet 

CX  Array  (50)  of  profile  velocitie*,  feet  per  second 

IPRNT  Debugging  print  flag 

COMMON  INPUTS 


/INIT/ 

I  ML 

Index  of  mixing  layer  in  profile  YX,CX 

PARAMETER 

OUTPUTS 

Yl 

Corrected  source  depth,  feet 

Y2 

Corrected  receiver  depth,  feet 

XRC 

Flag 

indicating  arrival  angle  depth 

COMMON  OUTPUTS 

/VELOC/ 

NPTSP 

No.  of  points  in  temporary 

YPP 

Array  (60)  of  corrected  profile  depths,  feet 

CPP 

Array  (60)  of  corrected  profile  velocities. 

feet  per  second 

/GRADS/ 

G 

Array  (60)  of  profile  gradients,  (mc"') 

/INIT/ 

K 1 

Index  of  ray  tracing  source  in  corrected 

prof  lie 

K2 

Index  of  ray  tracing  receiver  in  corrected 

profile 

YML 

Mixing  layer  depth  in  corrected  profile 

IMLP 

Index  of  mixing  layer  in  corrected  profile 

Cl 

Gradient  above  mixing  layer  in  corrected 

profile 

G2 

Gradient  below  mixing  layer  in  corrected 

profile 
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INSERT  (Cont'd) 


/CUSPCM/ 

ICUSP 

Cuaped  caustic  flag 

/PERIOD/ 

PERO 

Period  of  zero  degree  traced  from  Kl 

IAXFG 

Flag  indicating  aource-recciver-axia 

geometry 

COSA  ] 

^  Limiting  ray  angle  (at  axis)  for  analytic 

ANGP  J 

low  frequency  cutoff,  and  ita  coaine 

Begin  INSERT 


INSERT  1 


Continue  INSERT 


Is  source  depth  at  profile 
axis? 

'  No 


Are  source  and 
receiver  on  same 
side  of  profile  axis?j 
- [Yes  - 


Yes 


Bet  IAXFG-1  feet  IAXFG-2 


feet  IAXFG-1 


Call  RANGER  to  compute! 
period  of  2ero  degree 
ray  from  source 


If  mixing 
determine 
corrected  i 

Layer  exists 
Its  index  in 
irot i  le 

End  INSERT 


| 
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INSERT  2 


Yes 


Where  is  source? 

TZ - 

Between 


Tt - 

Surface 


5n 

Bottom 


Yes 


Is  source  depth  on  profile?} 

nm 


No 


Is  gradient 

In  isovelocity 

below  »  0? 
-  T  - - "J 

layer? 

Increase 

Increase 

velocity  at 

velocity  at 

profile  point 

end  of  layer 

by  2  feet/ 

by  2  feet/ 

second 

second 

No 


Yes 


Is  gradient 

Yes 

— « - ■ 

Is  gradient 

Yes 

below  •  0? 

above  -  0? 

(No 

f — - * - - - 

r 

Increase  source} 
depth  by  min. 
of  : 

a)  1/2  layer 
thickness 

b)  Distance  to 
get  2  ft/socj 
v@ locity 

c)  10  feet 


Decrease  source} 
depth  by  min. 
of : 

a)  1/2  layer 
thickness 

b)  Distance  to 
get  2  ft/sec} 
velocity 

:)  10  feet 
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INSERT  3 


SUBROUTINE  AXIS 


AXIS  is  called  from  INSERT  to  determine  the  parameters 
required  to  handle  axis-to-axis  transmission.  Initially, 
the  deep  sound  channel  axis  is  located  (if  present)  and 
the  smaller  of  the  two  velocities  immediately  to  either 
side  of  the  axis  is  obtained.  This  information  is  used 
to  fit  a  smooth  (quadratic)  function  to  the  velocity  pro¬ 
file  which  is  continuous  in  velocity  and  gradient  at  the 
axis.  The  period  of  the  2ero-degree  ray  at  the  axis  is 
then  determined  as  an  analytic  function  of  the  smoothed- 
profile  coefficients.  Subsequently,  the  angle  of  the  ray 
which  has  the  same  period  in  the  linearly-segmented  profile 
is  determined  by  means  of  >  sirple  closed-form  expression. 
If  the  two  turning-point  (vertex)  depths  of  this  ray 
about  the  axis  bound  both  the  source  and  receiver  depths, 
then  the  source  and  receiver  depths  are  required  to  take 
on  the  depth  value  of  one  of  the  turning  points. 

When  no  axis  exists,  the  “axis"  depth  index  is  set 
either  to  tho  surfaco  or  to  the  bottom,  and  the  period  of 
the  zero  degree  ray  (PERO)  Is  sot  to  0. 

Following  the  above  calculations,  AXIS  determines  the 
limit  angle  for  analytical  low-frequency  cut-off  effects. 
This  angle  is  the  angle  of  the  ray  at  the  axis  which  just 
grazes  the  next  layer  in  the  velocity  profile. 
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AXIS  (Cont'd) 


PARAMETER  INPUTS 


YS 

Source  depth,  feet 

YR 

Receiver 

depth,  feet 

IPR 

Debugging  print  flag 

COMMON  INPUTS 

/VELOCX/ 

NPTS 

No.  of  points 

in  sound  velocity  profile 

Y 

Array  (60)  of 

profile  depths,  feet 

C 

Array  (60)  of 

profile  velocities,  ft/sec 

PARAMETER 

OUTPUTS 

YS 

Adjusted 

source  depth, 

feet 

YR 

Adjusted 

receiver  depth 

i ,  f  eet 

COMMON  OUTPUTS 

/PERIOD/ 

YUP* 

Minimum  depth 

of  axis  ray,  feet 

YDN* 

Maximum  depth 

of  axis  ray,  feet 

PERO 

Period  of  axis 

ray,  feet 

COSA  ) 

Limit  angle  at 

axis  for  analytical 

ANGP  | 

lov-  frequency 

cutoff,  radians,  and  its  cosine 

LX 

Index  of  axis 

in  profile 

•LX  f  0  or  NPTS  only 
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SUBROUTINE  TABTH2 


TABTH 2  is  called  by  FACTTL  to  tabulate  the  ray  angle  at 
the  bottom,  as  a  function  of  the  ray  angle  at  the  receiver 
depth,  ©K2*  Tventy-one  equally-spaced  values  of  8^  are 
tabulated,  fro®  to  TT/2,  along  with  the  corresponding 

(unequally-spaced)  values  of  0_.  The  value  of  0...  is  that  of 
the  first  ray  which  touches  the  botto®.  For  ease  in  interpolation 
(by  function  THBOT} ,  the  ratios  of  the  corresponding  increments 
m  the  two  tables  are  also  tabulated. 

PARAMETER  INPUT 

CBC2  Ratio  of  botto®  velocity  to  receiver  depth 
velocity 


COMMON  OUTPUTS 

/TH2TAB/ 

TH2MIN 

Minimum  ray  angle  at  receiver  depth 

FACTOR 

Reciprocal  of  increment  in  8^ 

TH2T 

Array  of  21  values  of  8^ 

THBT 

Array  of  21  values  of  8 

D 

RATIO 

Array  of  20  ratios  of  A  0  /A  e 

B  3K2 

SUBROUTINE  CHITA 


CRITA  is  called  by  FACTTL  to  compute  the  KKB  phase 
factors  for  low-frequency  cutoff  effects.  The  deep-sound 
channel  is  first  located  (as  in  subroutine  AXIS),  and  sub¬ 
routine  RAYT  is  called  to  compute  the  relative  phase  of  the 
ray  along  the  axis  with  initial  angle  such  that  the  next  layer 
in  the  profile  is  just  grazed  (this  is  the  limit  angle  for 
analytical  cut-off  effects  which  was  calculated  by  AXIS) 

CRITA  then  determines  the  relative  phases  for  each  frequency 
being  processed,  using  analytical  expressions  for  rays  at  less 
than  the  limiting  angle,  and  by  iteration  (using  RAYT)  for  rays 
crossing  more  than  one  layer  in  the  profile.  For  each  frequency 
attenuation  factors  are  calculated  in  the  form  of  beam  patterns 
(amplitude  vs.  angle);  these  are  analytic  for  rays  below  the 
limit  angle  and  tabulated  f<->r  rays  above  the  limit  angle.  An 
array  of  flags  is  also  produced  for  rapid  determination  (in 
INSTOR  and  CUSP)  vt  which  beam  pattorn  type  applies  at  each 
f requency . 

PARAMETER  INPUT 

I PR  Debugging  print  flag 

COMMON  INPUTS 

/VELOC/  LHAX  No.  of  points  in  profile 

2  Array  (60)  of  profile  depths,  feet 

V  Array  (60)  of  profile  velocities,  ft/sec 
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CRITA  (Cont'd) 


/INIT/ 

L 

Index  of  source  depth  in  profile 

/RANGEL/ 

NFREQ 

No.  of  frequencies 

IFQMIN 

Index  of  lowest  frequency 

FQ 

Array  {6)  of  frequencies,  Hz 

/PERIOD/ 

COSA  1 

Limit  angle  for  analytical  beam  patterns, 

ANGP  ( 

radians,  and  its  cosine 

COMMON  OUTPUTS 

/CRIT/ 

BEE2 

Coefficient  of  analytic  low  frequency 

cut-off  amplitudes 

Cl 

Velocity  a-  source  ( K l > ,  ft/sec 

cx 

Velocity  at  profile  axis,  ft/sec 

JALF 

All  frequencies  -  analytical  flag 

JAIF 

Array  (6)  of  individual  analytical  frequency 

flags 

CRITANX 

Array  (6)  of  critical  angles  vs.  frequency 

CAX 

Array  (6,4)  of  beam  pattern  angles  vs. 

frequency 

SS 

Array  (6,4)  of  beam  pattern  amplitudes  vs. 

frequency 


SUBROUTINE  RAYT 


RAYT  is  called  irom  CRITA  to  determine  the  relative 
phase  (in  terms  of  travel  time)  along  a  ray  traced  over  one 
cycle,  adjusted,  if  required,  to  account  for  surface  and/or 
bottom  reflections. 

PARAMETER  INPUTS 

LA  Index  of  source  in  velocity  profile 

COSTHO  Cosine  of  initial  ray  angle 

COMMON  INPUTS 

/VLLOC/  LMAX  No.  of  points  in  velocity  profile 

Z  Arra>  (60)  of  profile  depths,  feet 

C  Array  (60)  of  profile  velocities,  ft/sec 

/GRADS/  G  Array  (60)  of  profile  gradients,  (sec  *) 

'ARAMETER  OCT  PIT 

CUTOFF  Relative  phase  of  ray,  seconds 


SUBROUTINE  ANGSCH 


ANGSCH  is  called  from  FACTTL  to  determine  the  families  of 
ray  i  to  be  processed,  as  a  function  of  the  sound  velocity 
profile,  and  the  source  and  receiver  depths.  Each  family  is 
chosen  so  that  a  smooth  fit  of  ray  range  vs.  ray  angle  can  be 
made.  Up  to  100  rays  in  up  to  20  families  are  allowed;  if 
profile  and  source/receiver  data  cause  these  maxima  to  be 
exceeded,  a  diagnostic  message  is  printed  (on  FORTRAN  unit  6) 
and  a  flag  is  set  to  indicate  this  condition.  The  ray  angles, 
m  radians,  are  stored  in  a  single  array,  THETA.  The  array 
IGRP  is  used  to  designate  the  index  of  the  first  ray  in  each 
family,  and  a  second  array,  1GRA2E ,  indicates  that  the  first 
ray  in  family  just  grazes  a  specified  profile  point. 

The  rays  m  each  group  are  constrained  to  be  at  leist 
three  in  number,  with  a  maximum  spacing  of  .5  degrees.  A  new 
family  begins  when  a  relative  maximum  in  the  profile  is 
encountered,  when  either  the  surface  or  bottom  is  encountered, 
and  when  the  profile  gradient  decreases  by  more  than  a  specified 
increment . 

COMMON  INPUTS 

/VELOC/  NPTS  No.  of  points  in  sound  velocity  profile 
Y  Array  (60)  of  profile  depths,  feet 

C  Array  (60)  of  profile  velocities,  ft/sec 
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ANGSCH  (Cont  *d) 


/GRADS/  G 

/IN IT/  Kl 

K2 
YKL 
IMLP 
I  ML 

/RAYZER/  THBINC 
THETMB 


Array  (60)  of  profile  gradient*, (sec**1 ) 

Index  of  source  depth  in  profile 
Index  of  receiver  depth  in  profile 
Mixing  layer  depth,  feet 
Index  of  mixing  layer  depth  in  profile 
Index  of  mixing  layer  (as  input) 

Bottom  angle  increment  (“5°),  radians 
Critical  angle  for  FNWC  bottom  type,  radians 


COMMON  OUTPUTS 
/RAYZER/  IRFFRZ 
/C'JSPCM/  ICUSP 

/ANGLES/  NRAYS 
NGRPS 
ISURF 
IBOT 
IGRP 
IGRAZE 
THETA 


Flag  indicating  grazing  ray  at  K2 
Flag  indicating  cusped  caustic  (turned  off 
if  source  and  receiver  within  mixed  layer) 
Number  of  rays  in  all  families  £  100  or  999 
Number  of  families  of  rays  £20 
Index  of  first  ray  of  hit  surface 
Index  of  first  ray  to  hit  bottom 

Array  (20)  of  indices  of  first  ray  in  each  family 

Array  (20)  of  grazing  flags 

Array  (100)  of  initial  ray  angles,  radians 
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SUBROUTINE  RANGER 


RANGER  is  called  from  FACTTL,  INSERT,  and  AXIS  to 
compute  the  ranges  associated  with  one  up-  and  down-going 
cycle  of  a  ray  traced  from  the  source  depth  (Kl) .  The  initial 
ray  angle  is  positive  if  possible  {source  not  at  the  surface) 
and  is  traced  for  a  quarter-cycle,  or  until  the  ray  hits  the 
surface;  the  process  is  then  repeated  for  a  ray  with  negative 
angle  if  possible  {source  not  at  the  bottom).  As  these  rays 
are  traced,  RANGER  calculates  depth  versus  range;  the 
crossings,  if  any,  of  the  receiver  depth  (K2)  are  noted, 
and  the  corresponding  ranges  are  saved. 


The  outputs  from  RANGER  are  the  period  of  the  ray,  the 
maximum  and  minimum  depths  attained  by  the  ray,  and  the  distances 
to  the  first  and  second  crossir.-js  s.  f  the  receiver  depth.  If 
the  receiver  depth  is  not  reached,  these  are  instead  the  ranges 
to  the  first  and  second  crossings  of  the  source  level,  i.e., 
the  half-  and  full-period  ranges. 


COMMON  INPUTS 
/VELOC/  NPTS 

y 

c 

/GRADS/  G 

/IN IT/  Kl 


No.  of  points  in  sound  velocity  profile 
Array  (60)  of  profile  depths,  feet 
Array  (60)  of  profile  velocities,  ft/sec 
Array  (60)  of  profile  gradients,  (sec  ') 
Index  of  source  depth 
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RANGER  (Cont'd) 


/ltilT/  K2  Index  of  receiver  depth 

AK  Vertex  velocity  of  ray  being  traced,  ft/sec 

SINTHO  Sine  of  initial  ray  angle  at  K1 

PARAMETER  OUTPUTS 

R1  Range  of  first  crossing  of  receiver  depth, 
f  eet 

R2  Range  of  second  crossing  of  receiver  depth, 
f  eet 

P  Full -cycle  period  of  ray,  feet 


COMMON  Ol'TPl’TS 

/PERIOD/  Yl’P  Minimum  depth  attained  by  ray,  feet 

YDN  Maximum  depth  attained  by  ray,  feet 
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SUBROUTINE  FITBOT 


FITBOT  is  called  by  FACTTL  to  determine  the  coefficients 

of  the  fit  of  range  vs.  ray  angle  at  the  receiver  (eK2)  *or  the 

last  (bottom-  and  surface-reflecting)  family.  For  rays  with 

o 

bottom  angles  less  than  30  ,  the  form  of  the  fit  is: 

R  -  Ain  .  a (2 ) . (®K3-®Hin> **  ♦  *(j).(eK2-eMln) 

For  bottom  angles  greater  than  30°,  the  form  of  the  fit  is: 

R  *  XTir*TAS7eK2 T TTT5T 

All  coefficients  are  calculated  by  FITBOT  with  the  exception  of 
A (4)  which  is  calculated  in  FACTTL  to  give  the  correct  limiting 
range  for  the  (implicit)  90°  ray.  The  derivatives  of  R  vs.  ©K2 
are  continuous  at  the  value  of  ©K2  corresponding  to  the  30° 
bottom  ray. 

PARAMETER  INPUTS 

R  Array  (2)  of  ranges  to  be  fit,  feet 

?H  Array  (2)  of  angles  (at  K2)  to  be  fit,  radians 

A  Array  value  (A (4)1  of  coefficient  determined  in  FACTTL 

PARAMETER  OUTPUTS 

A  Array  values  (a(1),  A(2),  A(3),  A(S))  of  coefficients 
of  fit  of  range  vs.  ray  angle 
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SUBROUTINE  FINDFT 


FINDFT  ia  called  by  FACTTL  to  determine  the  coefficients 
of  the  fit  of  ray  range  va.  ray  angle  at  the  receiver  for  all 
families  except  the  last  (bottom-  and  surface-reflected) 
family.  The  form  of  the  fit  is: 

R  -  A  ( 1 )  +  A(2UX(0k2)  ♦  A  ( 3)  •  (X(0K2))2 
The  form  of  the  function  X(0R2)  depends  upon  the  value  of 

I  PAR : 


I  PAR  ■ 

1: 

x(eK2> 

-  TAN  (6. 

I  PAR  - 

2: 

x(ex2) 

-  (eK2 

Three  points  are  used  in  the  fit:  The  minimum  and  maximum 
values  of  0R2>  plus  a  third  point  which  is  the  value  of  0R2  which 
gives  a  minimum  or  maximum  value  of  range  within  the  region  of 
interest,  or,  if  no  minimum  or  maximum  exists,  the  value  of 
0K2  corresponding  to  the  second  lay  in  the  family.  For  families 
with  grazing  rays  at  the  receiver  (Flag  IREFR2  »  1),  the  three 
values  of  0R2  are  0,  and  the  values  i@K2  corresponding  to  the 
last  ray  in  the  family. 

FINDFT  also  sets  various  parameters  giving  ranges  in 
range,  angle,  and  X,  over  which  the  fit  of  range  vs.  angle  is 
valid.  See  the  documentation  of  INSTOR  for  a  diagram  showing 
the  notation  employed  (page  4-83). 
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FINDFT  (Cont  *d) 


PARAMETER  INPUT 


IPRNT 

Debugging  print  flag 

COMMON  INPUTS 

/FLAGS/ 

C1C2 

Ratio  of  source  depth  velocity  to  receiver 

depth  velocity 

/INPUTS/ 

KRAY 

No.  of  rays  in  family  being  fit,  4  100 

KP 

Index  of  ray  paths  in  array  R,  1-4 

NORDEH 

Arrival  order  being  processed 

N  G 

Index  of  family  being  fit,  £  20 

I  PAR 

Flag  indicating  functional  form  of  fit 

R 

Array  '100,4)  of  ray  ranges  vs,  source 

angle  and  path  index,  feet 

TH 

Array  (100)  of  ray  source  angles,  radians 

/IN IT/ 

Kl 

Index  of  source  depth 

K2 

Index  of  receiver  depth 

/RAYZER/ 

IREFR2 

Flag  indicating  first  ray  grazes  receiver 

depth 

/CUSPCM/ 

I  CUSP 

Flag  indicating  cusped  caustic 

COMMON  OUTPUTS 

/flags/ 

THMIN 

Angle  at  receiver  of  shallowest  (at  source) 

ray,  radians 

THXAX 

Angle  at  receiver  of  steepest  (at  source) 

ray,  radians 

/FITS/ 

THF 

Array  (3)  of  receiver  ray  angles,  radians 
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FINDFT  (Cont ' d) 


/FITS/  RF  Array  (3)  of  ranges  of  fit,  feet 

A  Array  (3)  of  coefficients  of  fit 

XMIN  Value  of  X(THMIN) 

XMAX  Value  of  X(TKMAX) 

RMIN  Range  corresponding  to  THMIN 

RMAX  Range  corresponding  to  THMAX 

RANMIN  Minimum  range  spanned  by  fit,  feet 


SUBROUTINE  FITQ 


PITQ  is  called  by  FINDPT  and  CUSP  to  determine  the 
coefficients  of  the  quadratic  function 

Y  =  A  (1 )  ♦  A  (2 )  *X  ♦  A(3)*X2 

which  passes  through  the  three  points  X{1),Y(1);  X(2),Y{2); 
and  X (3) , Y (3) ,  where  X ( 1 )  *  X (2 )  +  X(3)  f  X(l). 

PARAMETER  INPUTS 

X  Array  (3)  of  ordinate  values 

Y  Array  (3)  of  abscissa  values 

PARAMETER  OUTPUTS 

A  Array  (3)  of  coefficients  of  fit 
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SUBROUTINE  INSTOR 


INSTOR  is  called  by  FACTTL  to  add  the  intensity  con¬ 
tributions  arising  from  one  arrival  of  one  order  of  one 
family  of  rays  to  each  applicable  point  in  array  TL  over 
the  range  of  interest.  INSTOR  is  also  called  by  CUSP  to 
process  any  smooth  caustic  which  may  be  associated  with  a 
cusped  caustic. 

The  range  interval  is  determined  by  the  coefficients 
and  parameters  of  the  fit  of  range  vs.  {receiver)  ray 
angle,  calculated  by  FACTTL  or  CUSP  and  passed  through 
common  areas  /FITS/  and  /FLAGS/.  The  contribution  of  each 
arrival  (which  may  actually  ref*esent  several  arrivals)  at 
a  given  range  point  within  this  interval  is  subsequently 
added  to  the  array  TL  by: 

7L ( Range , Freq !  *  TL { Range , Freqj  ♦  FACT/XL(Freq) 

XI  is  the  reciprocal  ray  intensity  and  is  determined  from 
one  of  four  arrival  geometries;  FACT  is  a  modifying  factor 
Ufi  the  range  0-4)  which  accounts  for  the  effects  of  in-, 
semi-,  or  f ully-coherent  combination  of  multiply-combined 
arrivals  under  the  ICON  and  IFLAG  options.  XI  may  also 
be  modified  to  account  for  bottom-bounce  losses  and  low- 
frequency  cut-off  effects. 

The  four  arrival  geometries  and  the  corresponding 
reciprocal  intensity  factors  are  as  follows.  See  relevant 
geometry  diagrams  on  pages  4-S3,  96  and  97. 


INSTOR  (Cont'd) 


1)  Single-ray  (no  caustic)  arrival: 


Xl(freq)  •  XITl  -  RANGE (K) *SINl ( 1 ) *RP/C0S2 


where ; 


Using : 


RANGE  (K)  «*  Range  in  £eet,R 


SINl(l) 
COS  2 
RP 


sine 


K1 


cos0K2i  (0K2 


(@K^  *  ray  angle  at  source) 

*  ray  angle  at  receiver) 


i 

JdR/d0K2j  at  range  R 


where : 


2)  Shadow  rone  of  a  snoot h  caustic: 

Xllfreq)  -  XIC ( f  r eg) *FA1RY  (XA1R) • RANGE (K;  /RC 

XXC(freq)  -  XICl*FREQ(freg) 

XIC1  -  CONST «S INC 1*RC* (ABS (RPP**2/SINC) ) •• 

(1 . / 3 . ) /C OSC 

XAIR  »  -BETA (f req) •ORAIR*FACTOR 

BETA ( f r  eq)  *  BETA 1 /FREQ (freq) **2 
DRAIR  -  ASS (RC -RANGE (K) ) 

FACTOR  -  1 . / (COS ( 1 . 570?9*AB$ (ORAIR/ (RCUT-RC) ) ) ) , 
*  !•  . 
BE7A1  -  (2, • (SINC/SGOO. ) *«2/A8$ (RPP) ) •* (l*/3. ) 

RANGE (K)  *  Range  in  feet,  R 

RC  •  Range  of  caustic,  feet 


FREQ 


(Radian  frequency) 


-1/3 


I RFC* 1 
IRFG-0 
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4)  Illuminated  region  of  a  smooth  caustic  associated  with 
a  cusp: 

a)  Steep  Branch  (j  0*  )  : 

Xl(freq)  is  calculated  a*  in  3)  above 

b)  Shallow  Branch  (J  eR2|<|  @c|  )  : 

Xl(freq)  -  AMINl (XI (freq) ,  XICP) 
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INSTOR  (Cont *  d ) 


Using : 

Xl(freq)  -  XIC(frcq)*FAIRY(XAIR)*2. * RANGE (K) /RC 
MCP  -  XI CUSP  (RANGE  (K)  ,1000.  *FREQK  ( f  req)  )  • RANGE  (K)  /COS2 

XlC(freq)  is  calculated  as  in  3)  above 

Where : 

FREQK(freq)  -  Frequency  in  Hr 

XI CUSP ( X)  •  Intensity  Function  -  see  CUSP,  page  4-93 

In  addition  to  calculating  the  intensity  contributions  as 
o,* lined  above,  INSTOR  determines  when  these  are  no  longer 
significant,  thus  signaling  the  end  of  the  loop  on  NORDEP  in 
FACTTL.  This  condition  arises  when  the  intensities  fall  below 
a  minimum  value,  or  when  the  range  of  an  order  is  greater  than 
the  maximum  range  of  interest.  Flag  IGTYP  in  common  area  /FLAGS/ 
is  set  negative  to  indicate  that  this  has  occurred. 

Arrival  information  (range,  ray  angle,  and  ray  intensities 
at  each  frequency)  is  calculated  and  written  on  file  IARVTP  if 
this  flag  is  not  zero. 

PARAMETER  INPUTS 

KRC  Flag  indicating  arrival  angle  depth 

SNTHR  Sign  of  ray  angle  at  depth  flagged  by  KRC 
IP  Debugging  print  flag 

IARVTP  Flag  Indicating  file  for  arrival  information  output 


INSTOR  (Cont  *d) 


PARAMETER  OUTPUTS 

TL  Array  (250,6)  of  intensities  vs.  range  and  frequency 
FILE  OUTPUT 

Unit  IARVTP  One  record  for  each  arrival  angle  at  each 

range  point.  The  format  of  each  BCD  record 
is  as  follows: 


Position 

Format 

Contents 

1 

1H 

♦ 

Blank 

2-7 

F6. 1 

Range,  Nautical  Miles 

8-9 

2X 

Blank 

10-16 

F7 . 3 

Arrival  Angle,  Degrees 

17-18 

2X 

Blank 

19-23 

F  5  .'-l 

Loss  at  1st  Frequency, 

dB 

24-25 

2X 

Blank 

26-30 

F  5 .  1 

Less  at  2nd  Frequency, 

dB 

etc 

etc 

etc 

52-58 

F5.  1 

Loss  at  6th  Frequency, 

dB 

COMMON  INPUTS 

/FLAGS/  IGTYPP 

Type  of  family  (fit) 

being  processed 

THMIN 

Shallowest  ray  angle 

in  family  being  processed 

THMAX 

Steepest  ray  angle  in 

family  being  processed 

CONST 

Coefficient  for  caustic  intensities 
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INSTOR  (Cont'd) 


C1C2  Ratio  of  source  velocity/receiver  velocity 

CUC2  Ratio  of  bottom  velocity/receiver  velocity 

ICOH  Flay  indicating  combination  of  arrivals 

IRSR  Flay  indicating  surface-reflected  rays 

SR 07  No.  of  bottom  bounces  of  family 

IRTU'  FNWC  bottom  type 

I  FLAG  Array  (6)  of  coherency  flags  vs.  frequency 

:k  Array  (2,6)  of  semi -coherent  phase  factors  vs 


f  r equenc y 

FITS/  A  Array  (5)  of  coefficients  of  fit  of  R  vs.  9 

XMIN  Minimum  value  of  argument  of  fit 

XMAX  M-ximum  value  of  argument  of  fit 

RMJN  Range  of  minimum-angle  ray  in  family 

RMAX  Range  of  maximum-angle  ray  in  family 

RANMIN  Minimum  range  at  which  intensities  result 
KP  Index  of  range  of  path  being  processed 

NCR PER  Order  of  arrival  of  path  being  processed 

N RANGE  No.  of  range  points 

NT  RFC  No.  of  frequencies 

JFQMIN  Index  of  minimum  frequency 

FNI u  Array  (6)  of  (radian  frequencies} **(- 1/3 ) 

FREQK  Array  (6)  of  frequencies,  KHz 

RANGE  Array  (250}  of  range  points,  feet 

/FAYZKP/  IRFFR2  Flag  indicating  grazing  arrival  at  receiver 


f  i  U  4  O 


/RANGEL. 
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INSTOR  (Confd) 


/AUTCOH/ 


/'CUSPCM/ 

/CRIT 


/PERIOD/ 


FNM1N 

FNMAX 

FNX I 
FNCYC 


ICUSP 

BEE2 

Cl 

CX 

JALF 

JAIF 

CRITANX 

CAX 

SS 

PERO 

ANGP 

IAXFG 


Min.  no.  of  range  points  per  surface- 
image  cycle 

Max.  no.  of  range  points  per  surface- 
image  cycle 

Reciprocal  of  FNMAX- FNMIN 
Array  (2.6)  of  cycles  of  phase  difference 
of  up-  and  down-going  rays  at  Kl  and 
K2  vs.  freq. 

Flag  ■  1  if  processing  smooth  caustic 
associated  with  a  cusp 
Coefficient  of  analytic  low-frequency 
cut-off  amplitudes 
Velocity  at  source  (Kl) ,  ft/sec 
Velocity  at  profile  axis,  ft/sec 
All-frequencies-analyt ical  flag 
Array  (6)  of  individual  analytical  frequency 
flags 

Array  (6)  of  critical  angles  vs.  frequency 
Array  (6.4)  of  beam  pattern  angles  vs. 
frequency 

Array  (6,4)  of  beam  pattern  amplitudes  vs. 
frequency 

Period  of  zero-degree  ray  at  source,  feet 
Limit  angle  for  analytical  beam  patterns 
Flag  indicating  source-receiver-axis  geometry 
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INSTOR  (Cont'd) 


COMMON  OUTPUTS 

/FITS/'  RANMIN  Minimum  range  at  which  fit  is  applicable 

FLAGS/  lGTYPi  Set  negative  to  indicate  no  contribution 

to  intensity 


4-82 


R  is  quadratic  in  X1 (0)  or  X2 (©) 
i.e.,  R  -  A  ( 1 )  4  A(2)*X<0)  +  A(3)*[X{0))2 


Xx (0)  -  TAN  (0)  IPAR  -  1 

x,(6)  *  e  -  «  IPAR  -  2 

*  Min 


INSTOR  -Fit  of  R  vs.  6 

NOTE:  Throughout  this  figure  e  •  eR2  «  Ray  angle  at  receiver  (K2) 
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4  SNTHP  ! 

J 


I  2  +1 .  -1.  |  - 

\  ,  i 

a  * 

I  2  SNTHR  SNTHP 

- - < - 1 - ■ - 

4  :  -1.  i  -1.  i  *1. 

.  L. _ _  _ , _ 1  - . . i _ .  j 


Determination  of  No.  of 
Arrival  Angles  and  Their  Signs 
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RCtrr  -  [factor  in  table  (n-norder)J  *pero 


IREFR2  -  1  so  KP  3,4 


ICL'SP  “  0: 
IAXFG  -  1: 
KP  - 


ICAUST 


IAXFG  -  2: 


KP  « 


INSTOR  -  Determination  of  RCUT 


Begin  INSTOR 


INSTOR  1 


Continue  IN5TOR 


4-87 


Continue  INSTOR 


INSTOR  3 
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Continue  INSTOR 


Exit  INSTOR 
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INSTOR  i 


It 


In  shadow  zone  of  caustic? 


In  illuminated  region  of 
caustic? 

Ives  f  No 


Beyond  maximum  range 
of  expression  of 
caustic  intensity? 

JVes 

^No 

Processing 

Processing 

shal low 

shallow 

brarch  from 

branch 

CU3p  • 

arising 

[no 

from  cusp? 

j  Tfes 

1 

r 

Yes  i  No 

Calculate  in¬ 
tensity  from 
one- ray 
expression 


Calculate  in¬ 
tensity  from 
Ilium,  caustic 
expression 


TTTculate  In¬ 
tensity  from 
maximum  of 
Ilium,  zone 
and  cusp  in¬ 
tensity 


Calculate  in¬ 
tensity  from 
shadow  sore 
caustic  ex¬ 
pression 


Attenuate  intensity  if 
ray  has  undergone  bottos 
bounces 


INSTOR  $ 
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:  V- 


Continue  INSTOR  -  Loop  on  frequency 


Determine  factor  for  modify- 
inq  intensity  to  account  for 
coherency  geometry  and  for 
coherency  addition  type 


Modify  intensity  by  coherency 
and  low-frequency  cutoff 
factors 


Add  intensity  to  TL  array  as 
a  function  of  range  and  fre¬ 
quency 
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SUBROUTINE  CUSP 


CUSP  is  called  by  FACTTL  to  add  the  intensity  contributions 
arising  from  each  order  (other  than  direct)  of  a  family  containing 
a  cusped  caustic,  to  each  applicable  point  in  the  array  TL  over 
the  range  oi  interest.  CUSP  also  calls  INSTOR  to  process  any 
smooth  caustic  which  may  be  associated  with  a  cusped  caustic. 

The  range  interval  is  determined  directly  from  ray  path 
distances  for  the  cusp  itself;  CUSP  also  calculates  the 
coefficients  and  parameters  of  the  fit  of  range  vr.  (receiver) 
ray  angle  for  the  smooth  caustic  (if  any)  which  are  passed 
through  common  area  /FITS/  to  INSTOR.  For  each  range  point  in 
the  interval  covered,  the  intensity  contribution  to  each  frequency 
arising  from  one  arrival  order  is  added  to  the  array  TL  by: 

TL(Range,  F req)  ~  TL(Range,  Freq)  ♦  1. /XI? (Freq) 

XIP  is  the  reciprocal  ray  intensity  and  may  be  modified  to  account 
for  bottom-bounce  losses  and  iow-f requeney  cutoff  effects,  it 
is  computed  from: 

XIP(Treq)  «  XI (Freq) *RAN/CST 

Using : 

XI(Freq)  •  XICUSP (IR, RAN, 1000. sPREQIMFreq) ) 

CST  -  COS  (THETA) 

Where : 

RAN  -  Range  in  feet,  R 
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CUSP  (Cont'd) 


THETA  ■  Angle  of  ray  at  range  R,  estimated  from 
linear  or  quadratic  fit  depending  upon 
ray  being  processed. 

XICUSP  *  Function  computing  vertical  spreading  for 

a  ray  near  the  cusp  as  a  function  of  region 
(one  or  three  ray) ,  range,  and  frequency. 

IR  ■  1  in  one-ray  region,  ?  in  three-ray  region 

of  cusp. 

In  additior  to  calculating  cusp  intensity  contributions  as 
outlined  above,  CUSP  determines  when  these  are  no  longer  significant, 
thus  signaling  the  end  of  the  loop  on  NORDER  in  FACTTL.  This 
condition  arises  when  the  intensities  from  the  smooth  cauptics 
calculated  by  INSTOR  fall  below  a  minimum  value,  or  when  the 
range  of  an  ordor  is  greater  than  the  maximum  range  cf  interest. 

p  i  t  fvM1 17  i  e  —  ..  j  :  i  —  —  j  j  i.  i  , 

&  *yv« c#  tT*  s%0\  tu  imutait;  uub  vynuu  iu,i. 

Arrival  information  (range,  ray  angle,  and  ray  intensities 
at  each  frequency)  is  calculated  and  written  on  file  IARVTP  if 
this  flag  is  not  zero. 

PARAMETER  INPUTS 

Kl  Index  of  source  depth 

X2  Index  of  receiver  depth 

KRC  Flag  indicating  arrival  angle  depth 


CUSP  (Cont 'd) 


IPRNT  Debugging  print  flag 

IARVTP  Flag  indicating  file  for  arrival  information 

output 


PARAMETER  OUTPUTS 

I DONE  Flag  indicating  end  of  contribution  to  intensities 
f L  Array  (250,6)  of  intensities  vs.  range  and 

frequency 

FILE  OUTPUT 

Unit  1ARVTP  Or.e  record  for  each  arrival  angle  at  each  range 
point.  See  INSTOR  for  format  specifications. 


COMMON  INPUTS 


/INPUTS/ 

N* 

R 

TH 

/RANGEL/ 

N RANGE 

NFREQ 

FREQK 

RANGE 

/FLAGS/ 

IGTYP 

NBOT 

IBTYP 

/CUSPCM/ 

CCUSP 

/GRADS/ 

G 

Number  of  rays  in  family 
Array  (100,4)  of  ranges  vs.  ray  angle  and 
arrival,  feet 

Array  (100)  of  ray  angles  of  family,  radians 

No.  of  range  points 

No.  of  frequencies 

Array  (6)  of  frequencies ,KHs 

Array  (250)  of  ranges,  feet 

Type  of  family  being  processed 

No.  of  bounces  at  this  order 

FNWC  bottom  type 

Velocity  at  cusp,  ft/sec 

Array  (60)  of  profile  gradients*  ;«0C*’t) 


•KRAYS  elsewhere 
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CUSP  (Cont'M! 


/PERIOD/ 

AN’GP 

/CHIT/ 

BEE  2 

Cl 

cx 

JALF 

JAIF 

CRITANX 

CAX 

SS 

COMMON  OUTPUTS 
/FITS/  THF 

RF 

A 

XMIN 

XMAX 

RMIN 

RMAX 

RANMIN 

/C US PCM/  RCUSP 

BCUSP 


Limit  angle  for  low-frequency  cutoff 
Coefficient  of  low-frequency  cut-off  amplitudes 
Velocity  at  source  (Kl) ,  ft/sec 
Velocity  at  profile  axis,  ft/«ec 
All-f requenciea-analytical  flag 
Array  (6)  of  individual  analytical  frequency 
flags 

Array  (6)  of  critical  angles  vs.  frequency 
Array  (6,4)  of  beam  pattern  angles  vs. 
f  requency 

Array  (6,4)  of  beam  pattern  amplitudes  vs. 
frequency 


Array  (3)  of  angles  of  fit  oi  R  vs.  0  for 
smooth  caustic,  radians 
Array  (3)  of  ranges  of  fit,  feet 
Array  (3)  of  coefficients  of  fit 
Minimum  value  of  argument  of  fit 
Maximum  value  of  argument  of  fit 
Range  of  minimum-angle  rsy  in  family 
Ranee  of  maximum-angle  ray  in  family 
Minimum  range  at  which  intensities  result 
Range  of  cusp,  feet 
Cusp  parameter  beta 
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Ray  Region  Three-Ray  Region 


CUSP  -  Parameters  Of  Fit  Of  Smooth  Caustic 


CUSP  1 


Continue  CUSP 


Loop  over  range  point*:  K-KMIN,  UMAX 


i 


Determine  no.  of  rays  at 
point  «*  NRAYS 

range 

1 

Compute  intensity  factor 
frequency 

for  each 

1 

Loop  over  number  of  rays: 

1-1,  NRAYS 

i 


Determine  angle  of  ray 
being  processed  and  corres- 
ponding  cut-off  angle _ 


1 


Loop  over  frequency: 

NF-1,  NFREQ 

_J 

Compute  intensity 
contribution  of  ray, 
attenuate  if  ray  has 
undergone  bottom 
bouncos,  modify  by 
low-freq.  cutoff 
factor 

Add  intensity  to  TL 
array  as  a  function  of 
range  and  frequency 

L _ _ _ j 

_ _ 1 

Output  arrival  angle  inform*-* 
tion  if  flagged 


A 


Continue  CUSP 


— 

Smooth  caustic  also  present? 

Yes 

Exit 

No 

CUSP 

. 

Determine  coef'is.  of  fit 
type  1  and  associated 
parameters 

J 

Call  SETSNR  to  determine 
arrival  angle  signs 

L 

Call  INSTOR  to  process  smooth 
caustic 

.  ..  I  . .  ... 

Did  intensity  ot  min-ranga 
smooth  caustic  fail  to 
contribute  significantly 

Yes 

'  T 

Set  flag  to  indicate  in¬ 
tensity  contribution  made 

Exit 

> 

CUSP 

J _ 

Set  flag  to  indicate  no 
intensity  contribution 

Exit 

CUSP 

CUSP 
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SUBROUTINE  HFCHTL 


HFCHTL  is  called  by  FACTTL  to  calculate  the  intensities 
for  other  than  direct  and  bottom-surface-ref lected  paths  when  the 
half-channel  indicator  has  been  aet  on  input  (index  of  mixing 
layer  at  bottom  depth:  IML  ■  NPTS) .  The  routine  ia  atrictly 
applicable  only  for  three  source/receiver  depth  combinationa 
(60/60  feet,  200/300  feet,  200/60  feet)  and  frequencies  (50, 

300,  850,  1700  Hz).  The  computed  intenaitiea  are  proportional 
to  1/Range;  the  constant  of  proportionality  ia  initially 
expressed  as  a  transmission  loss  (dB  re  1  foot)  in  the  form: 

TL  -  A  -  B  •  log  (D/1000) 
where  D  is  the  ha  1 f -channel  depth  in  feet. 

The  coefficients  A  and  B  are  each  chosen  from  a  separate 
3x4  array  as  a  function  of  source/receiver  depth  combination  and 
of  frequency.  The  coefficients  in  these  arrays  were  themselves 
computed  by  FACT,  using  a  temporary  correction  aet  to  bypass  the 
contributions  to  transmission  loss  arising  from  the  direct  and 
bottom-surface-ref lect ing  paths.  These  FACT  runs  were  made 
using  a  simple  pressure-gradient  profile  Uvelocity/Adepth  • 

-  i 

.018  sec  )  to  half-channel  depths  of  1000  and  18,000  feet. 

PARAMETER  INPUTS 

NR  No.  of  range  points,  £  250 
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HFCHTL  (Cont'd) 


NPREQ 

No.  of  frequencies,  £  6 

VS 

Source  depth,  feet 

YR 

Receiver  depth,  foet 

D 

Dopth  of  half-channel,  feet 

FREQ8 

Array  (6)  of  frequencies,  Hr 

RANGE 

Array  (250)  of  range  points, 

PARAMETER  OUTPUTS 

TL  Array  {2  SO, 6)  to  which  intensities  are  added 

as  a  function  of  range  and  frequency. 
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SUBROUTINE  QUAD 

QUAD  is  called  from  INSTOR  to  solve  the  quadratic  equation 
expressing  R  as  a  function  of  X(0),  in  order  to  find  the  value(s) 
cf  X (8)  (and  thus  the  value(s)  of  8),  if  any,  at  the  range  point 
being  processed  in  INSTOR.  The  roots  are  constrained  to  lie 
between  XMIN  *  x^eMin^  and  XMAX  *  X^Max*‘  r°turns  the  number 

of  roots  *  0,  1  or  2,  and  the  corresponding  (ordered)  values  of 
X (8) .  See  diagram  page  4-81. 

The  notation  in  QUAD 

Y  *  A(l)  +  A  ( 2 ) *X  +  A  (3)  *X**2 

corresponds  to 

R  -  A(l)  +  A  ( 2 )  •  X  ( 0 )  4-  A  (  3)  *X  (0)  **2 

in  INSTOR. 

PARAMETER  INPUTS 

A  Array  (3)  of  coefficients  of  fit 
XMIN  Minimum  value  allowable  for  X 
XMAX  Maximum  value  allowable  for  X 
Y  The  range  for  which  X  values  are  desired 

PARAMETER  OUTPUTS 

X  Array  (2)  of  solutions  of  equation 
IFLAO  No.  of  X  values  returned  (0,  1  or  2) 
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FUNCTION  SPEED 

SPEED  ia  called  from  INSERT  to  linearly  interpolate  in  the 
sound  velocity  profile  to  determine  the  velocity  corresponding  to 
a  source  or  receiver  depth  which  has  not  been  explicitly  included 
as  a  point  on  the  input  profile.  The  results  are  unpredictable 
if  this  depth  is  less  than  the  minimum  profile  depth;  the  last 
segment  of  the  profile  is  linearly  extrapolated  for  a  depth 
greater  than  the  maximum  profile  depth. 

PARAMETER  TNPUT 

VP  Depth  at  which  velocity  is  to  bo  determined,  feet 

COMMON  INPUT 

/VELOCX/  NPTS  No.  of  points  on  profile 

Y  Array  (60)  of  profile  depths,  feet 

C  Array  (60)  of  profile  velocities,  ft/sec 

FUNCTION  OUTPUT 

SPEED  Interpolated  value  of  velocity  at  depth  YP,  ft/sec 
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FUNCTION  SET8NR 


SETSNR  is  called  from  FACTTL  and  CUSP  to  determine  the  sign<s) 
of  the  'n  ival  angle  (a)  at  the  depth  at  which  arrival  infor¬ 
mation  is  being  determined.  It  return*  the  sign(s)  of  the 
angle  of  the  ray  (a)  corresponding  to  the  path  of  interest#  which 
may  be  multiple,  according  to  the  coherency  option  in  effect. 

If  KRC  -  1,  the  e-ign  of  the  angle (s)  at  depth  K1  is  returned; 

if  KRC  -  2,  the  sign  of  the  angle (a)  at  depth  K2  is  returned. 

A  value  of  0.  is  returned  to  indicate  that  both  up-  and  down¬ 
going  rays  (signs  ■  +1.  and  -1.)  are  present. 

PARAMETER  INPUTS 

ICOH  Coherency  option  in  effect  -  0,1,2  or  3 

KRC  Flag  indicating  arrival  depth 

Ki  Index  of  source  depth 

K2  Index  of  receiver  depth 

KP  Index  of  path  being  processed 

FUNCTION  OUTPUT 

SETSNR  SignU)  of  arrival  angles:  ♦!.,  -1.,  0. 

See  below  for  values  returned  ss  a  function  of 
input  parameters. 
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FUNCTION  FAIRY 


FAIRY  is  called  by  INSTOR  s*  a  step  in  calculating  tha 
intensities  which  exist  in  the  neighborhood  of  a  smooth 
caustic.  The  function  is  related  to  the  Airy  function  (for 
real  arguments)  as  follows: 


FAIRY  IX) 


-a. 77  i-  -X  ^4.0 


X  C  -4.0 


-  .7*2  •  (X)i/2  X  >  1.77 


The  function  is  approsimatad  over  the  primary  range  of 
interest  by  10.**CX,  where  CX  is  interpolated  from  tabulated 
values  of  C ( I )  vs.  I;  these  valuas  are: 


C  (I ) 


2  log 


10 


Ai  (0) 


Ai(^^-) 


X  •  1(1)31 


This  range  of  I  corresponds  to  -X  ranging  from  -2.  to 
♦4.  at  intervals  of  0.2. 


PARAMETER  INPUT 

X  Argument  of  function 

FUNCTION  OUTPUT 


FAIRY  Functional  value  corresponding  to  X 
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FUNCTION  XICUSP  _ 

XICUSP  ie  called  by  CUSP  and  INSTOR  aa  a  step  in 
calculating  the  intenaitiea  which  exiat  in  the  neighborhood 
of  a  cusped  cauatic.  The  function  computes  the  vertical 
spreading  near  the  cusp,  returning  the  aingle-path  value 
assuming  three  equal-amplitude  paths  in  the  interference 
(three-ray)  region  and  one  path  in  the  one-ray  region. 

The  calculation  is  throe-step: 


1)  X  - 


T4^ 

_  c  J 


n  1/2 


where 


m  Cusp  parameter 
W  °  Radian  frequency 
C  ■  Cusp  sound  velocity 


2)  y  -  f[*(R-Rc)  •  X) 

whore : 

R  -  Range  of  interest,  feet 
r^  «  Range  of  cusp,  feet 
f  is  calculated  by  call  to  function  FE2B.  The 
sign  of  the  argument  is  positive  in  the  one- 
ray  region,  negative  in  the  three- ray  region 


3)  XICUSP  •  •  V) 


PARAMETER  INPUTS 

IR  Pleg  indicating  one-  or  three-ray  region 
ran  Range  of  interest,  feet 


FREQ  Frequency,  H* 
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XICUSP  (Cont'd) 


FUNCTION  OUTPUT 


XICUSP 


Functional  value  of  epreading 


FUNCTION  PE2B 

PE2B  is  called  from  XICUSP  to  calculate  the  Pearcey- 
function  component  of  the  expression  for  vertical  spreading 
in  the  vicinity  of  a  cusp.  The  function  is  related  to  the 


Pearcey  function  as  follows: 

P£2B(y)  -*ir/-y 

y  -i.  5 

-  Pe  (y) 

-3.5  <y  <0 

-  Pe(y) 

0.  <  y  <2. 

-  TP/y 

y  >  2.0 

The  function  Pe(y)  is  approximated  over  the  primary  rang? 
of  interest  by  P**2,  where  P  is  interpolated  from  tabulated 
values  of  PT(I)  vs.  I;  I  ■  1(1)23.  This  range  of  I  corresponds 
to  y  ranging  from  -3.5  to  2.0  at  intervals  of  0.25. 

PARAMETER  INPUT 

Y  Argument  of  function 

FUNCTION  OUTPUT 

PE2B  Functional  value  corresponding  to  y 
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FUNCTION  THBOT 

THBOT  is  called  by  INSTOR  and  CUSP  to  determine  the  ray 
angle  at  the  bottom  as  a  function  of  the  ray  angle  at  the 
receiver  depth.  This  angle  is  determined  from  linear  inter¬ 
polation  in  the  values  of  bottom  angle  vs.  receiver  ray  angle 
as  previously  tabulated  by  routine  TABTHZ. 


PARAMETER  INPUTS 

TH2  Ray  angle  at  receiver  depth  (K2) 


COMMON  INPUTS 
/TH2TA8/  TH2MIN 

FACTOR 
TH2T 
THAT 
RATIO 


Minimum  ray  angle  at  receiver  depth 
Reciprocal  of  Increment  of  TH2T  values 
Array  (21)  of  ray  angles  at  receiver  depth 
Array  (21)  of  ray  angles  at  the  bottom 
Array  (20)  of  ratio  of  increments  of 
TH2T  to  increment  of  THBT 


FUNCTION  OUTPUT 

TttBQT  Interpolated  value  of  bottom  angle 
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FUNCTION  BOTTOM 


BOTTOM  is  called  from  INSTOR  and  CUSP  to  calculate  the 
intensity  attenuation  occurring  along  bottom-bounce  ray  paths. 

The  attenuation  is  returned  as  a  value  between  0  and  1. 

Currently,  BOTTOM  returns  a  1  {no  attenuation)  ir  the  input 
bottom  type  is  0,  and  calls  upon  FNVfC  routine  BTMLOS  for  bottom 
types  1-9.  This  routine  may  be  replaced  or  restructured  according 
to  the  desires  of  the  FACT  user. 

PARAMETER  INPUTS 

NBOT  No.  of  bottom  bounces 

1BTYP  Bottom  type,  0-9 

THB  Ray  angle  at  bottom,  radians 

FR£0  Frequency,  Ht 

FUNCTION  OUTPUTS 

BOTTOM  Bottom  attenuation 
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5. 


THE  FACT  HANDOUT 


The  following  pages  contain  a  reproduction  of  the  AESD 
FACT  Handout,  a  computer-maintained  document  which  is  an 
integral  part  of  the  FACT  Package.  The  various  sections  are 
as  follows: 

•  A  description  of  the  physical  basis  of  the  FACT  Model} 
e  A  description  of  the  computer  program  flowi 
e  A  description  of  the  FACT  Package  Program  Library; 
e  A  description  of  the  input  card  formats  accepted  by 
the  TLOSS  program; 

e  A  description  of  the  input  deck  structure; 
e  A  listing  of  a  sample  set  of  input  data  cords; 
e  The  outputs  resulting  from  this  input  deck; 
e  The  listing  produced  by  a  CDC  6600  FORTRAN  compilation 
Of  FACT. 

Of  the  last  section,  only  the  listing  of  the  program  TLOSS 
is  included  in  this  report. 
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DEPARTMENT  OF  THE  NAVY 

OFFICE  OF  NAVAL  RESEARCH 
875  NORTH  RANDOLPH  STREET 
SUITE  1425 

ARLINGTON  VA  22203-1995 


IN  REPLY  REFER  TO: 


5510/1 

Ser  321OA/01 1/06 
31  Jan  06 


MEMORANDUM  FOR  DISTRIBUTION  LIST 

Subj :  DECLASSIFICATION  OF  LONG  RANGE  ACOUSTIC  PROPAGATION  PROJECT 
(LRAPP)  DOCUMENTS 

Ref:  (a)  SECNAVINST  5510.36 

Enel:  (1)  List  of  DECLASSIFIED  LRAPP  Documents 

1 .  In  accordance  with  reference  (a),  a  declassification  review  has  been  conducted  on  a 
number  of  classified  LRAPP  documents. 

2.  The  LRAPP  documents  listed  in  enclosure  (1)  have  been  downgraded  to 
UNCLASSIFIED  and  have  been  approved  for  public  release.  These  documents  should 
be  remarked  as  follows: 

Classification  changed  to  UNCLASSIFIED  by  authority  of  the  Chief  of  Naval 
Operations  (N772)  letter  N772A/6U875630,  20  January  2006. 

DISTRIBUTION  STATEMENT  A:  Approved  for  Public  Release;  Distribution  is 
unlimited. 

3.  Questions  may  be  directed  to  the  undersigned  on  (703)  696-4619,  DSN  426-4619. 


BRIAN  LINK 
By  direction 


Subj :  DECLASSIFICATION  OF  LONG  RANGE  ACOUSTIC  PROPAGATION  PROJECT 
(LRAPP)  DOCUMENTS 

DISTRIBUTION  LIST: 

NAVOCEANO  (Code  N121LC  -  Jaime  Ratliff) 

NRL  Washington  (Code  5596.3  -  Mary  Templeman) 

PEO  LMW  Det  San  Diego  (PMS  181) 

DTIC-OCQ  (Larry  Downing) 

ARL,  U  of  Texas 

Blue  Sea  Corporation  (Dr. Roy  Gaul) 

ONR  32B  (CAPT  Paul  Stewart) 

ONR  3210A  (Dr.  Ellen  Livingston) 

APL,  U  of  Washington 

APL,  Johns  Hopkins  University 

ARL,  Penn  State  University 

MPL  of  Scripps  Institution  of  Oceanography 

WHOI 

NAVSEA 

NAVAIR 

NUWC 

SAIC 


Declassified  LRAPP  Documents 


